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Insights into the nature of northwest-to-southeast aligned 
ionospheric wavefronts from contemporaneous Very Large Array 
and ionosondes observations 

J. F. Helmboldt 1 



The results of contemporaneous summer nighttime obser- 
vations of midlatitude medium scale traveling ionospheric 
disturbances (MSTIDs) with the Very Large Array (VLA) 
in New Mexico and nearby ionosondes in Texas and Col- 
orado are presented. Using 132, 20-minute observations, 
several instances of MSTIDs were detected, all having wave- 
fronts aligned northwest to southeast and mostly propagat- 
ing toward the southwest, consistent with previous studies 
of MSTIDs. However, some were also found to move toward 
the northeast. It was found that both classes of MSTIDs 
were only found when sporadic- F (F s ) layers of moderate 
peak density (1.5 < foEs < 3 MHz) were present. Limited 
fbEs data from one ionosonde suggests that there was a 
significant amount of structure with the E s layers during 
observations when foEs> 3 MHz that was not present when 
1.5 < foEs < 3 MHz. No MSTIDs were observed either be- 
fore midnight or when the F-region height was increasing at 
a relatively high rate, even when these E s layers were ob- 
served. Combining this result with AE indices which were 
relatively high at the time (an average of about 300 nT 
and maximum of nearly 700 nT), it is inferred that both 
the lack of MSTIDs and the increase in F-region height are 
due to substorm-induced electric fields. The northeastward- 
directed MSTIDs were strongest post-midnight during times 
when the F-region was observed to be collapsing relatively 
quickly. This implies that these two occurrences are related 
and likely both caused by rare shifts in F-region neutral 
wind direction from southwest to northwest. 



1. Introduction 

Traveling ionospheric disturbances (TIDs) are one of the 
most well known/studied ionospheric phenomena. So-called 
medium-scale TIDs, or "MSTIDs," are some of the most 
commonly observed TIDs with wavelengths and speeds on 
the order of 100 km and 100 m s _1 , respectively. They have 
been thoroughly studied, mainly with GPS-based data [e.g., 
Hernandez- Paj ares et a/., 2006; Tsugawa et a/., 2007] and 
airglow imagers [e.g., Shiokawa et a/., 2003, 2008], primarily 
at midlatitudes where they are among the most common im- 
pact ive disturbances. They are particularly problematic for 
high-precision GPS-based applications requiring centimeter- 
level accuracy [Hernandez- Paj ares et al., 2006]. 

It has been demonstrated that during nighttime, espe- 
cially within the summer, the wavefronts of MSTIDs in the 
northern hemisphere are almost exclusively aligned north- 
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west to southeast, typically propagating toward the south- 
west [Hernandez- Paj ares et al, 2006; Tsugawa et a/., 2007]. 
The origin of this orientation and direction is still not firmly 
established. In the past, the orientation of these structures 
was attributed to the Perkins instability Perkins [1973]. 
However, the direction of propagation toward the southwest 
in this case is the opposite of what one would naively ex- 
pect given the typical eastward-directed background ExB 
drift. This was addressed by Kelley and Makela [2001] who 
showed that with a northeastward-directed Pedersen cur- 
rent present resulting from a southeastward-directed neutral 
wind, the finite extent of the MSTID wavefronts leads to a 
polarized electric field pointing along the wavefronts toward 
the northwest. Thus, ExB drift from this polarized field 
will push the wavefronts toward the southwest. 

Despite the qualitative agreement with observations, the 
Kelley and Makela [2001] model suffers from the drawback 
that it relies on the Perkins instability, which has a growth 
rate that is too low to generate these features [~ 10 -4 
s _1 ; Cosgrove et al, 2004]. By demonstrating that the 
Perkins instability is coupled to a separate instability within 
sporadic- F (F s ) layers, Tsunoda and Cosgrove [2001] and 
Cosgrove et al. [2004] showed that northwest-to-southeast 
aligned features can form within both the F-region and 
the E s layer with a combined growth rate that is roughly 
2-5 times larger. Since the formation of these structures 
within the F-region is still driven, at least in part, by the 
Perkins instability, the same polarized electric fields pre- 
dicted by Kelley and Makela [2001] should form along the 
wave fronts. For the optimum configuration, Cosgrove et 
al. [2004] showed that the horizontal velocity of structures 
within the F-region is proportional to — F n e + F e n, imply- 
ing that the northwestward directed, polarized electric field 
predicted by Kelley and Makela [2001] would cause south- 
westward motion. 

While this theoretical framework provides a plausible de- 
scription of how nighttime, midlatitude MSTIDs are gen- 
erated, there remains little direct evidence to support it. 
In addition, three-dimensional simulations of coupled dy- 
namics between E s layers and the F-region showed that in 
some cases, polarized electric fields form perpendicular to 
the northwest-to-southeast aligned wavefronts rather than 
parallel [Yokoyama et al, 2009]. In these cases, southwest- 
ward motion of the phase fronts was attributed to the south- 
ward neutral wind with the coupled E s layer which was also 
demonstrated by the simulations of Cosgrove [2007]. The 
Cosgrove [2007] simulations also showed that the strength of 
the coupling was reduced when a significant amount of struc- 
ture was present within the E s layer. Recent studies based 
on observations of disturbances over Japan [Shiokawa et al., 
2008] and the southwestern United States [Helmboldt et al., 
2012a] have also shown that summer nighttime MSTIDs can 
propagate toward the northeast, sometimes reversing direc- 
tion from previously being directed southwest ward. These 
northeastward-directed waves have been tentatively linked 
to significant drops in F-region height, which may provide 
a clue to the mechanisms that determine the direction of 
motion of nighttime, midlatitude MSTIDs. 
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To gain further insights into the nature of these distur- 
bances, contemporaneous data from the Very Large Ar- 
ray (VLA) in New Mexico (specifically, 34° 04' 43.497" N, 
107° 37' 05. 8 19" W) and nearby ionosondes have been com- 
piled, spanning 17 nights in June 2002. Using techniques pi- 
oneered by Helmboldt and Interna [2012], VLA observations 
of multiple cosmic sources were used to characterize the en- 
vironment of ionospheric fluctuations while using ionosodes 
data to constrain the properties of both E s layers and the F- 
region. The results are presented in §2 and discussed further 
in §3. 

2. Fluctuation Spectra, Sporadic- and F- 
region Compression 

To explore the relationship between MSTIDs and both 
E s layers and F-region height, a data-set of contemporane- 
ous VLA and ionosondes observations during summer night- 
time was compiled. The VLA data were taken from a larger 
database generated by the VLA Low-frequency Sky Survey 
[VLSS; Cohen et a/., 2007], a 74 MHz survey of the north- 
ern sky conducted between 2001 and 2007. The VLSS con- 
sisted of three or more 20-minute observations for each of 
523 overlapping, circular fields, each 15° across, yielding a 
total of 1862 separate observations. The observations were 
conducted with the VLA interferometer in its B configura- 
tion, an inverted "Y" of 27 dishes spanning approximately 
11 km. The correlated signals from pairs of antennas, or 
"visibilities," were used to image each 15° field. During the 
imaging process, the affects of the ionosphere across each 
field were compensated for as a function of time using 1-2 
minute "snapshot" images of between 20 and 50 relatively 
bright cosmic sources. A process referred to as field-based 
calibration [Cotton et a/., 2004] uses the position shifts of 
these sources to map the effects of the ionosphere over the 
entire field. 

Each source position shift is proportional to the mean 
gradient in the total electron content (TEC) over the span 
of the array along the line of sight to the source. Given the 
typical rms noise within a single snapshot, the TEC gradi- 
ents measured from these position shifts have a precision of 
about 10~ 4 TECU km" 1 [1 TECU= 10 16 m" 2 ; Helmboldt 
and Interna, 2012]. Thus, these gradient measurements can 
be used to probe a much larger range in TEC fluctuation am- 
plitudes than is possible with similar GPS-base data which 
yield a typical precision on the order of 0.01 TECU [e.g., 
Hernandez- Paj ares et al., 2006]. 

It was shown by Helmboldt and Interna [2012] that the 
mean TEC gradient time series derived from cosmic source 
position shifts can be used with a straightforward Fourier 
analysis to produce a three-dimensional spectral cube of 
TEC gradient fluctuations, one temporal dimension, and 
two spatial. This analysis consists of first computing a sep- 
arate power spectrum cube for each component of the gra- 
dient (i.e., north-south and east- west) using temporal and 
spatial discrete Fourier transforms. For each spectrum, the 
impulse response function (IRF) is computed using the po- 
sitions of the observed sources and the spectrum is normal- 
ized by the amplitude of the IRF. This way, for a simple 
plane wave, the spectral power will be just the square of its 
amplitude. However, since not all observed structures are 
wavelike, the IRF normalization is denoted in the units re- 
ported for the spectra, typically (mTECU km -1 IRF -1 ) 2 , 
where mTEC is a "milli-TECU" (i.e., 10" 3 TECU). 

Next, the power spectra for the two gradient components 
are added together so that for a plane wave with spatial 
frequency £ and TEC amplitude A, the combined spectral 



power will just be (2tt^A) 2 . Note, however, that (1) the 
spectra are convolved with the IRF which has a full width 
at half maximum of about 0.009 km -1 and (2) many of the 
observed fluctuations are not wavelike. Therefore, the com- 
bined spectra do not necessarily go to zero at the origin. For 
instance, turbulent fluctuations have a power-law spectrum 
approximately proportional to £ -11 ^ 3 for TEC fluctuations 
and £ _5//3 for the combined TEC gradient spectra [see, e.g., 
Helmboldt et a/., 2012a]. This spectral analysis has been ap- 
plied to all 1862 VLSS observations, and a complete clima- 
tological investigation based on the resulting spectral cubes 
is presented by Helmboldt et al. [2012b]. 

All of the summer observations within the VLSS were 
conducted in June 2002 (under VLA program number 
AP441). The mean nighttime (local times between 20:00 
and 04:00) spectral cube from these observations presented 
by Helmboldt et al. [2012b] shows a prominent detection of 
southwestward- directed MSTIDs, as expected. For refer- 
ence, the mean spectral power over all temporal frequencies 
from this average spectral cube is displayed in Fig. 1 as a 
function of north-south and east-west spatial frequencies, 
^ns and ^ew- Since there were 132 summer nighttime ob- 
servations, this provides an excellent opportunity to inves- 
tigate how the spectral properties of northwest-to-southeast 
aligned TEC fluctuation wavefronts depend on other iono- 
spheric properties. In particular, given the discussion in §1, 
it is interesting to investigate how the spectra change with 
E s properties and with the height of the F-region. 

To this end, data from nearby ionosondes that were ob- 
tained during the June 2002 VLSS observations were com- 
piled. The two ionosondes closest to the VLSS observations 
were at Boulder, Colorado and Dyess Air Force Base (AFB) 
in Texas. The locations of these stations relative to the iono- 
spheric pierce-points for the 132 summer nighttime VLSS 
observations, assuming a nominal ionospheric height of 300 
km, are shown in Fig. 2. One can see that the Boulder 
station is at a similar longitude, but higher latitude that 
the VLSS observations; the opposite is true for the Dyess 
AFB station. Together the two stations provide a reason- 
able depiction of ionospheric conditions near the VLSS ob- 
servations. 

To constrain the relevant E s and F-region properties, val- 
ues for the maximum E s reflection frequency, foEs, and the 
F-region height, h'F, were obtained from both stations for 
each VLSS observation. Unfortunately, the E s blanketing 
frequency, fbEs, was not available from the Dyess AFB sta- 
tion and was flagged as unreliable most of the time (67 out of 
132 observations) for the Boulder station during the VLSS 
observations. Therefore, the VLSS observations could not 
be binned according to either fbEs or foES— feEs without 
losing more than half of the observations. However, the 
limited Boulder fbEs measurements can be useful for inter- 
preting the results presented below and will be revisited in 
§3. 

For each VLSS observation, the closest ionosondes ob- 
servations within ±20 minutes were used. For foEs, 84 out 
of 132 VLSS observations had measured values from both 
ionosondes stations, for which the mean foEs between the 
two stations was used. The remaining 48 had data from 
Dyess AFB only. For h'F, several measurements from the 
Boulder station were flagged as unreliable and were not used, 
leaving 89 VLSS observations with h'F from Dyess AFB 
only. 

During all of the 132 summer nighttime VLSS observa- 
tions, E s was detected with at least one of the two ionoson- 
des stations. The measured foEs from each station ranged 
from about 1.5-6 MHz. However, there were several in- 
stances where E s was detected by one station and not the 
other. In these cases, foEs=0 was used for non-detections, 
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lowering the mean foEs below 1 MHz. The summer night- 
time VLSS observations were consequently put into four bins 
based on foEs ranging from 0-6 MHz and a mean fluctua- 
tion spectrum cube was computed for each bin. The mean 
spectral power over all temporal frequencies is displayed for 
each bin in the upper row of Fig. 3 as a function of ^ns and 
^ew- From these, one can see that both southwestward- 
directed and northeastward-directed MSTIDs are seen only 
in the 1.5<foEs<3 MHz bin. Note that the weaker exten- 
sion visible in the northeastern quadrant for this foEs bin is 
not prominent within the mean summer nighttime spectrum 
shown in Fig. 1. 

The results found by Shiokawa et al. [2008] and Helm- 
boldt et al. [2012a] discussed in §1 provide motivation to 
explore how the direction of propagation of these MSTIDs 
depends on the change in F-region height. This was done 
by first estimating the h'F temporal gradient, dh'F/dt, for 
each VLSS summer nighttime observation. This was done 
by numerically computing the temporal derivative using La- 
grangian interpolation for each of these observations with at 
least two other VLSS observations within ±40 minutes. This 
excluded only one of the 132 nighttime observations. It was 
found that most (54%) of the observations had dh'F/dt be- 
tween —30 and 30 km hr -1 . Roughly equal fractions had 
dh'F/dt < -30 km hr" 1 (20%) and dh'F/dt > 30 m hr" 1 
(26%). 

The observations with 1.5 < foEs < 3 MHz were split 
into three bins, one with dh'F/dt < —30 km hr -1 , one with 
-30 < dh'F/dt < 30 km hr" 1 and one with dh'F/dt > 30 
km hr -1 . A map of mean spectral power is displayed in 
the middle row of Fig. 3 for each of these bins. The re- 
sults are rather striking and show that for dh'F/dt < —30 
km hr -1 , the northeastward-directed waves are nearly as 
strong as those directed toward the southwest. Both classes 
of MSTIDs are weaker for the -30 < dh'F/dt < 30 km hr -1 
bin with the southwestward-directed waves as the dominant 
wavelike feature. Finally, neither class of waves in present 
within the spectral map for the <9h'F/<9t>30 km hr -1 bin. 

The southwest and northeast features in the dh'F/dt < 
—30 km hr -1 spectral map seem so symmetric that one may 
suspect they are the result of an artifact within the Fourier 
analysis. However, when examined along the temporal fre- 
quency axis, these two features are quite distinct, peaking 
at two different temporal frequencies. This is shown in Fig. 
4 where maps from the mean spectral cube for temporal 
frequencies of v — 7.5 hr -1 and v — 10.5 hr -1 are plotted 
which represent the peak frequencies for the southwest and 
northeast features, respectively. The data used to make the 
mean spectral cube for this 3\\F /dt bin were also carefully 
examined and no reason was found to doubt that these fea- 
tures are real detections of separate phenomena. 

The spectra from observations with 1.5 < foEs < 3 MHz 
were also binned by local time into pre-midnight and post- 
midnight bins. The mean spectral maps for these two bins 
are displayed in the bottom panels of Fig. 3. One can 
see that the pre-midnight spectrum is very similar to the 
dh'F/dt > 30 km hr -1 spectrum in that it shows no signs of 
wave activity. Both southwest and northeast features can be 
seen in the post-midnight map with the southwest feature 
being the dominant component. Several combinations of lo- 
cal time bins were tried for the post-midnight observations, 
but none produced spectral maps with northeastward fea- 
tures as strong as that seen in the map for dh'F/dt < —30 km 
hr -1 . This indicates that the northeastward-propagating 
MSTIDs are in fact preferentially seen during times when 
the F-region height is dropping relatively quickly. 

3. Discussion 

While not immediately obvious, the results presented in 
Fig. 3 are consistent with many aspects of the theoretical 
framework presented in §1 when other ionospheric condi- 
tions are taken into account. First, the fact that sum- 
mer nighttime MSTIDs were only seen when E s layers were 



present is qualitatively consistent with these waves being 
generated by the E-F coupling instability (EFCT) described 
by Cosgrove et al [2004]. 

That these waves were preferentially seen when moder- 
ate rather than high values of foEs where measured may 
indicate that E-F coupling is generally weaker when there 
is a larger degree of structure within the E s layer. This 
presupposes that the larger observed values of foEs origi- 
nate from dense structures within the E s layers rather than 
from uniform E s layers which are simply more dense. This 
is supported by the fbEs measurements available from the 
Boulder station for 65 of the 132 summer nighttime VLSS 
observations. Fig. 5 shows that among the Boulder sta- 
tion data, foEs is near or below the blanketing frequency 
for observations where foEs< 3 MHz during which all of the 
MSTID activity was found. In contrast, for most of the ob- 
servations where foEs> 3 MHz, the difference between foEs 
and fbEs exceeds 1 MHz. Ogawa et al. [2002] noted that 
radar echoes from ^-region structures on scales of tens of 
kilometers or less [i.e., quasiperiodic (QP) echoes] generally 
appear when foEs— fbEs> 1 MHz. They also showed that 
the strength of these QP echoes correlates with foEs— fbEs. 
Thus, when considered together, the results shown in Fig. 
3 and 5 strongly support the notion that an increased level 
of structure within E s layers compromises their ability to 
interact electrodynamicly with the Perkins instability in the 
F-region. This is consistent with results of the simulations of 
E-F coupling presented by Cosgrove [2007] that were men- 
tioned in §1. 

Among observations with moderate foEs values, observa- 
tions that are either in the pre-midnight sector or are as- 
sociated with larger than normal upward F-region motions 
show no signs of MSTIDs. This can be explained by taking 
general conditions during the time of the observations into 
account. During the 132 summer nighttime VLA observa- 
tions used, the mean AE index was 276 nT, indicating that 
substorm activity was likely common. In addition, the mean 
among pre-midnight observations was 315 nT with a maxi- 
mum of about 690 nT, as apposed to a mean and maximum 
of 198 nT and 390 nT, respectively, for their post-midnight 
counterparts. 

Park [1971] noted that during substorm activity, zonal 
electric fields are commonly induced within the ionosphere 
which typically switch from eastward pre-midnight to west- 
ward post-midnight. The collapse in F-region height caused 
by this shift is often seen near midnight as a negative tempo- 
ral gradient in TEC within VLA data colloquially referred 
to as the "midnight wedge" [see, e.g., Kassim et al., 2007; 
Helmboldt and Interna, 2012, Fig. 14 and 11, respectively] 
Thus, it is not surprising that the spectral map for obser- 
vations with large upward F-region motions is similar to 
those that were conducted before midnight since this is the 
most likely time for eastward electric fields to exist during 
substorm activity. In addition, Cosgrove et al. [2004] noted 
that if the background electric field is (magnetic) eastward 
with little or no north-south component, the F-region is sta- 
ble against the Perkins instability. Thus, one would expect 
to see little if any wave activity pre-midnight if the EFCl, 
which depends on the Perkins instability, is chiefly respon- 
sible for the growth of northwest-to-southeast aligned wave- 
fronts. In contrast, the "resonant" condition for the EFCl 
given by Cosgrove et al. [2004] has an electric field with 
a westward component. This would imply that wave activ- 
ity should appear post-midnight when the substorm-induced 
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electric field switches from eastward to westward, which is 
precisely what was observed. 

In the post-midnight sector during times when the F- 
region height is relatively stable, the observations show ev- 
idence of MSTIDs that are almost exclusively directed to- 
ward the southwest. At midlatitudes in June, Emmert et 
al. [2003] showed that the neutral winds in the F-region are 
usually directed toward the southwest, especially after mid- 
night. This implies that it is unlikely that the polarized elec- 
tric field described by Kelley and Makela [2001] contributes 
significantly to the southwestward motion of these waves 
since it requires a southeastward-directed neutral wind. It 
is more plausible that the background electric field, which 
was likely westward during the post-midnight VLSS observa- 
tions, was influencing the motion of these waves as predicted 
by Cosgrove et al. [2004]. However, the typical southward 
neutral wind velocity within the E s layers present at the 
time could also have had a significant effect, similar to what 
was seen in the simulations of Yokoyama et al. [2009] . 

The fact that northeastward-propagating waves were only 
seen when the F-region was dropping at a relatively high 
rate implies that these drops in height are related to the di- 
rection of motion for these waves. As suggested by Shiokawa 
et al. [2008], the northeastward-directed MSTIDs may be 
associated with rare instances when the neutral wind shifts 
direction toward the northwest. In this case, one would 
expect polarized electric fields to form along the MSTID 
wavefronts, similar to that described by Kelley and Makela 
[2001], but in the opposite direction. Thus, the motion in- 
duced by these fields via ExB drift would be toward the 
northeast. In addition, the northward component of the neu- 
tral wind would help the westward electric field drive down 
the F-region height at an even higher rate, which matches 
what was observed. 

Note that, in order for this scenario to work, the influence 
of the polarized electric fields along the wavefronts would 
have to be stronger than those of both the background west- 
ward electric field and the southward neutral wind within 
the E s layers. This may only be true during some fraction 
of the time when northwestward neutral winds are present. 
Therefore, even when the F-region is collapsing at a higher 
than normal rate, which may be the signature of a northward 
meridional wind, the wavefronts could propagate toward the 
southwest. The equally strong presence of southwestward- 
directed waves during times when dh'F/dt < —30 km hr _1 
seems to support this. Three-dimensional simulations sim- 
ilar to those performed by Yokoyama et al. [2009] with a 
northwestward neutral wind could potential provide further 
insights. 

Finally, note that the results presented here are hampered 
somewhat by the fact that the VLSS observations were con- 
fined to a single month and year during a period when sub- 
storm activity was apparently common. A survey of cosmic 
sources in the northern sky is being planned that will use 
the new 330 MHz system being developed for the VLA by 
the Naval Research Laboratory and the National Radio As- 
tronomy Observatory. Unlike the VLSS, this survey will be 
scheduled to optimize time-of-day and time-of-year cover- 
age. In doing so, it is also likely to sample a wider range 
in levels of geomagnetic activity than the VLSS observa- 
tions presented here. It will also benefit from contempora- 
neous GPS data from nearly 40 new, continuously operat- 
ing receiver stations throughout New Mexico that will pro- 
vide the ability to better characterize medium to large-scale 
ionospheric structures. In addition, a new digisonde station 
will be operating at Kirtland AFB in nearby Albuquerque. 
Merging these three data-sets will lead to the development 
of a comprehensive picture of midlatitude ionospheric dis- 
turbances, including E s and MSTIDs. 



Acknowledgments. Basic research in astronomy at the 
Naval Research Laboratory is supported by 6.1 base funding. The 
VLA is operated by the National Radio Astronomy Observatory 
which is a facility of the National Science Foundation operated 
under cooperative agreement by Associated Universities, Inc. The 
author is grateful to the UK Solar System Data Centre from which 
the ionosondes data was obtained and the World Data Center 
for Geomagnetism, Kyoto for providing AE indices. The author 
would also like to thank the reviewers for useful suggestions which 
greatly improved the paper. 



References 

Cohen, A. S., W. M. Lane, W. D. Cotton, N. E. Kassim, T. J. 
W. Lazio, R. A. Perley, J. J. Condon, and W. C. Erickson 

(2007) , The VLA Low- Frequency Sky Survey Astron. J., 134, 
1245-1262 

Cotton, W. D., J. J. Condon, R. A. Perley, N. Kassim, J. Lazio, 
A. Cohen, W. Lane, and W. C. Erickson (2004), Beyond the 
isoplanatic patch in the VLA Low-frequency Sky Survey, Proc. 
SPIE, 5489, 180 

Cosgrove, R. B., R. T. Tsunoda, S. Fukao, and M. Yamamoto 
(2004), Coupling of the Perkins instability and the sporadic E 
layer instability derived from physical arguments, J. of Geo- 
phys. Res., 109, A06301 

Cosgrove, R. B. (2007), Generation of mesoscale F layer struc- 
ture and electric fields by the combined Perkins and E s layer 
instabilities, in simulation, Ann. Geophys., 25, 1579-1601 

Emmert, J. T., B. G. Fejer, and D. P. Sipler (2003), Climatology 
and latitudinal gradients of quiet time thermospheric neutral 
winds over Millstone Hill from Faby-Perot interferometer mea- 
surements, J. Geophys. Res., 108 (A5), 1196 

Helmboldt, J. F. and H. T. Interna (2012), Very Large Array Ob- 
servations of Disturbed Ion Flow from the Plasmasphere to the 
Nighttime Ionosphere, Rad. Sci., in press 

Helmboldt, J. F., T. J. W. Lazio, H. Interna, and K. Dymond 
(2012), A New Technique for Spectral Analysis of Ionospheric 
TEC Fluctuations Observed with the Very Large Array VHF 
System: From QP Echoes to MSTIDs Rad. Sci., RS0L02 

Helmboldt, J. F., W. M. Lane, and W. D. Cotton (2012), Clima- 
tology of Mid-latitude Ionospheric Disturbances from the Very 
Large Array Low-frequency Sky Survey, Rad. Sci., submitted 

Hernandez-Parajes, M., J. M. Juan, and J. Sanz (2006), Medium- 
scale traveling ionospheric disturbances affecting GPS mea- 
surements: Spatial and temporal analysis, J. of Geophys. Res., 
Ill, A07S11 

Kassim, N. E. T. J. W. Lazio, W. C. Erickson, R. A. Perley, 
W. D. Cotton, E. W. Greisen, A. S. Cohen, B. Hicks, H. R. 
Schmitt, and D. Katz (2007), The 74 MHz System on the Very 
Large Array, Astrop. J. Supp., 172, 686-719 

Kelley, M. C, and J. J. Makela (2001), Resolution of the discrep- 
ancy between experiment and theory of midlatitude F-region 
structures, Geophys. Res. Let., 28, 2589-2592 

Ogawa, T., O. Takahashi, T. Otsuka, K. Nozaki, M. Yamamoto, 
and K. Kita (2002), Simultaneous middle and upper atmo- 
sphere radar and ionospheric sounder observations of midlati- 
tude E region irregularities and sporadic E layer, J. Geophys., 
107(A10), 1275 

Park, C. G. (1971), Westward Electric Fields as the Cause of 
Nighttime Enhancements in Electron Concentrations in Mid- 
latitude F Region, J. Geophys. Res., 76, 4560-4567 

Perkins, F. (1973), Spread F and Ionospheric Currents, J. Geo- 
phys. Res., 78, 218-226 

Shiokawa, K., Y. Otsuka, C. Ihara, T. Ogawa, and F. J. 
Rich (2003), Ground and satellite observations of nighttime 
medium-scale traveling ionospheric disturbances at midlati- 
tude, J. Geophys. Res., 108(A4), 1145 

Shiokawa, K., Y. Otsuka, N. Nishitani, T. Ogawa, T. Tsugawa, T. 
Maruyama, S. E. Smirnov, V. V. Bychkov, and B. M. Shevtsov 

(2008) , Northeastward motion of nighttime medium-scale trav- 
eling ionospheric disturbances at middle latitudes observed by 
an airglow imager J. of Geophys. Res., 113, A 123 12 



HELMBOLDT: VLA AND IONOSONDES OBSERVATIONS OF MSTIDS 



X- 5 



Tsugawa, T. Y., Y. Otsuka, A. J. Coster, and A. Saito (2007), 
Medium-scale traveling ionospheric disturbances detected with 
dense and wide TEC maps over North America, Geophys. Res. 
Lett, 34, L22101 

Tsunoda, R., R. B. Cosgrove (2001), Coupled electrodynamics in 
the nighttime midlatitude ionosphere, Gephys. Res. Let., 28, 
4171-4174 



Yokoyama, T., D. L. Hysell, Y. Otsuka, and M. Yamamoto (2009), 
Three-dimensional simulation of the coupled Perkins and Ens- 
laver instabilities in the nighttime midlatitude ionosphere J. 
of Geophys. Res., 114, A03308 



J. F. Helmboldt, US Naval Research Laboratory, Code 
7213, 4555 Overlook Ave. SW, Washington, DC 20375 
(joe. helmboldt @nrl . navy, mil) 



HELMBOLDT: VLA AND IONOSONDES OBSERVATIONS OF MSTIDS 



mean summer nighttime VLSS spectrum 



00 W 



-3.15 



-3.30 



-3.45 - 



-3.60 



-3.75 



0) 

o 

CL 



-3.90 



-4.05 f 



-4.20 



-0.02 



0.00 

e w km" 1 



0.02 



Figure 1. The mean TEC gradient fluctuation spectrum 
over all temporal frequencies for all summer nighttime 
(20:00-04:00 local time) VLSS observations. 
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Figure 2. The ionospheric pierce points for the VLSS 
observations conducted during nighttime (20:00-04:00 lo- 
cal time) in June 2002, assuming an ionospheric height of 
300 km (blue points). The locations of the two ionoson- 
des used are plotted as x 's. 
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Figure 3. Upper: The mean spectral power (over all 
temporal frequencies) from observations binned by foEs 
as a function of north-south and east-west spatial fre- 
quency, ^ns and £ew> For observations in the second 
foEs bin, 1.5<foEs<3 MHz, mean spectral power maps 
with the observations binned by h'F temporal gradient, 
dh'F/dt, and local time are displayed in the middle and 
lower rows of panels, respectively. 
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Figure 4. For observations where 1.5 < foEs < 3 MHz 
and dh'F/dt< —30 km hr -1 (see the lower middle panel 
of Fig. 3), maps of the spectral power at the peak tempo- 
ral frequencies of the southwestward (left; 7.5 hr -1 ) and 
northeastward (right; 10.5 hr -1 ) propagating waves. 
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Figure 5. For the 65 observations with available 
measurements for foEs and fbEs from the Boulder sta- 
tion (fbEs was not available from the Dyess station), 
foEs— fbEs versus foEs. The foEs bins used to make the 
mean fluctuation spectra shown in the upper panels of 
Fig. 3 (excluding the lowest foEs bin) are shaded sepa- 
rately for reference. 



